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ABSTRACT: The proof of the stability of matter is three decades old, but the 
question of stability when arbitrarily large magnetic fields are taken into account was 
settled only recently. Even more recent is the solution to the question of the stability of 
relativistic matter when the electron motion is governed by the Dirac operator (together 
with Dirac's prescription of filling the "negative energy sea"). When magnetic fields are 
included the question arises whether it is better to fill the negative energy sea of the free 
Dirac operator or of the Dirac operator with magnetic field. The answer is found to be 
that the former prescription is unstable while the latter is stable. 



I. Introduction 

This will be a brief summary of recent work which solves the problem of the stability 
of matter in magnetic fields. The first part, which was the one reported in the Evora 
conference, concerns the usual Schrodinger equation with nonrelativistic kinetic energy, 
but with a Zeeman term cr ■ B. The second part is about some very recent work [1], not 
reported at the conference, on the relativistic problem with kinetic energy given by the 
Dirac operator. In both cases the electrons are treated quantum mechanically while the 
fields are regarded as classical fields. 

It has been well understood for some time, starting with the work of Dyson and Lenard 
in 1967 [2], that matter consisting of electrons and K static (but arbitrarily positioned) 
nuclei, and governed by nonrelativistic quantum mechanics, is stable. (See also [3] and see 
[4], [5] for reviews.) This means that the ground state energy (more generally, the bottom 
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of the spectrum, in case there is no lowest eigenvalue) is finite and is bounded below by a 
universal constant times the number of particles, i.e., Eq > const. {N + K). The reason the 
nuclei are taken to be static is that they are so massive compared to the electron that if 
the nuclear mass (or the nuclear radius of 10~^^cm, for that matter) played any significant 
role then ordinary matter would have to look very different from what it does. 
The Hamiltonian of our system in suitable units is 

N 

-f^nonrelatvistic = ^ Pi + aVc , (1.1) 

where a = /%c is the fine-structure constant and p = — iV, as usual. is the Coulomb 
potential of K fixed nuclei with nuclear charge Ze (we take all the charges to be the same 
for simplicity only), with locations Rj in R'^, and with N electrons. 

N K 

= -Z^^\xi- Rj\-^ + \xi-Xj\-^ + Z^ \R,-Rj\-\ (1.2) 

i=l j=l l<i<j<N l<i<3<K 

We should like to generalize this stability question in two directions. One direction 
is to add magnetic fields and the other is to make the kinetic energy relativistic, which 
means replacing by ^/p"^ + m? — m or, equivalently and more simply by |p| (since the 
difference between these two expressions is bounded by m, and hence the change in energy 
is at most a term proportional to the number of electrons). It will turn out that these 
two problems (magnetic fields and relativity) are not as unrelated mathematically as one 
might think. Later on, we shall generalize still further by replacing the operator |p| by the 
Dirac operator; this is the content of the very recent work mentioned above. 

II. Stability of non-relativistic matter with magnetic 
fields 

Let us begin with the magnetic field problem, and with a simple example to illustrate 
its essence. We shall consider the nonrelativistic case first. We shall denote by A{x) 
the vector potential of a magnetic field and B(a;) = curl A(a;) is the magnetic field itself. 
Conceptually, we do not ask where this field comes from; it can be external or it can be 
generated by the current of the electrons themselves. The action of the field on the orbital 
motion of the electrons is given by the minimal substitution (recalling that the electron's 
charge is negative) 

p^p + A(a;). (2.1) 

This has little effect on the ground state energy of our Hamiltonian in (1.1). Indeed, for 
one electron the field can only raise the energy. For several electrons it can lower the 
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energy but by a controllable amount. In fact, all known proofs of the stability of matter 
hold with constants unchanged by the substitution (2.1). 

The problem comes in when we include the electron spin-field interaction. Thus, we 
replace (p + A)^ by the Pauli operator (in which a denotes the triplet of Pauli matrices) 

Ta = [a • (p + A)]2 = {p + Af + a-B. (2.2) 

It is a fact, which is not entirely easy to prove [6], that Eq for a hydrogen-like atom 
H = 7a — Z/r, behaves as — (InS)^ for a large, spatially uniform field B. Thus, if we 
minimize Eq over all possible fields (since nature will presumably find that value of the 
field that gives the lowest Eq) the minimum energy will be —oo. However, we forgot to 
include the self energy of the magnetic field 

^^fieid = [87ra]-^ J B{xf d^x. (2.3) 

(The unconventional factor in (2.3) comes from the system of units we have been 
using.) Does the addition of -fffleid stabilize the system, i.e., is Eq now bounded below, 
independent of B? 

The answer, perhaps not unsurprisingly in view of the appearance of a in (2.3), is 
yes if and only if a is small enough. More precisely it is the combination Za^ that has 
to be small enough. This was proved in [7], where it was shown that the "hydrogenic" 
Hamiltonian 

Hhydrogenic = ^A — ^o/r + -fffield (2.4) 

is bounded below if Za^ < Qtt^/S. Oddly, it was harder to prove that the atom is unstable 
if Za^ is too large; this fact relied on a difficult calculation in [8] that the Pauli kinetic 
energy operator 7a can have a zero mode for certain special choices of the A field. That 
is, the equation 

Ta* = (2.5) 

has a solution for some magnetic field B, whose energy i/fieid is finite, and some square 
integrable, electron wave function . The ordinary kinetic energy does not have this 
remarkable property, i.e., the property that in a suitable field the electron can "stand still". 
(This happens for the uniform magnetic field, to be sure, but that field has infinite energy.) 

Obviously, the next question is whether the full many-electron, many-nucleus problem 
is stable when is replaced by 7a and with the addition of i?fieid- Thus, we consider 

N 

i^Pauli = ^Ta {i) + OiV^ + Afield , (2.6) 
i=l 

and ask whether it is bounded below (for all A, AT, K and all locations of the nuclei) by 
some universal constant times N + K. This problem was open for about 9 years and was 
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first solved by Fefferman (unpubfished) for sufficiently small a, Z (which does not include 
hydrogen in its range). A simple proof of stability with good constants was given shortly 
thereafter in [9]: A representative sufficient condition for stability in (2.6), obtained in [9], 
is 



Another sufficient condition is Z < 1050 when a = 1/137. (A more recent result in [1] 
improves this to Z < 2265.) It is comforting to have this very large margin of stabil- 
ity in (2.6) because it means that the magnetic field-electron interaction can be treated 
perturbatively with confidence when Z < 92. 

It is important to understand why there is a bound on a alone in (2.7), independent of 
the magnitude of Z (as long as Z > 0). The reason is that once one admits the possibility 
that a single atom can collapse if Za^ is too large (or Za is too large in the relativistic 
theories studied below) then one has to ask what prevents a large number of small nuclei 
from coming together to form a supercritical nucleus? The answer is the Coulomb repulsion 
among the nuclei, but this quantity is effectively governed by the parameter (or in 
the relativistic theory) because the nuclear repulsion is proportional to Z^a — {Zc^Y 1^ ■ 
In other words if Zo? is truly the relevant atomic parameter then it should be regarded 
as fixed, which means that oT^ is the relevant parameter governing the nuclear repulsion. 
Anyway, whether or not the reader believes this heuristic discussion, the fact remains that 
the bound on ct in (2.7) is essential. 

Two proofs of stability are given in [9]. The simplest uses an ingredient that, at first 
sight appears to have nothing to do with the magnetic problem — the stability of relativistic 
matter. Let us turn to that next. 

III. Relativistic matter without magnetic fields and its 
application to the nonrelativistic magnetic problem 

The relativistic Hamiltonian looks like (1.1), with the substitution (2.1), but with 
replaced by |p|. It contains no i^fleid term. Thus, 



Admittedly, this Hamiltonian is not well motivated physically because the operator 
|p| is not local (unlike the Dirac operator). Minimal substitution does preserve gauge 
invariance, however. There is no spin-field interaction, as in the Pauli or Dirac Hamiltoni- 
ans, but let us take one step at a time. It will turn out that the understanding of (3.1) is 
central to the whole story. 

There are two important remarks to be made about (3.1). The first is that everything 
scales like 1/length (if we let the A field scale as 1/length as well). Because of this, there 



Za^ < 0.041 and a < 0.06. 



(2.7) 



N 




(3.1) 
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are exactly two possibilities for £'0 after we minimize over all possibilities for A and the 
nuclear locations: 

Eo = or Eo = -00 . (3.2) 

(If we used y^p"^ + — m instead, then the first alternative, Eq = would be replaced 
by some negative number depending on m. The condition on a for the first alternative to 
hold does not depend on m.) 

The second remark is that for the one-electron hydrogenic problem i?o = if and only 
if Za <2/tt (due to Kato [10] and Herbst [11]). The reason for such a condition is that the 
kinetic and potential energies scale the same way. If the kinetic energy is larger, Eq can 
be driven to zero by dilation; in the reverse case, Eq can be driven to —00 by contraction. 
Closely related to this is a property of the Dirac hydrogenic atom that was known for a 
long time: Self adjointness of the Dirac operator requires Za < 1. 

The stability of the many-body Hamiltonian (for all and K) was first shown by 
Conlon [12] for Z = 1 and a < 10"^'^'^ and zero magnetic field. The situation was sub- 
stantially improved in [13] to Za < 2/tt and a < 1/94, but A = 0. [13] also contains a 
simpler proof of stability which holds only up to a slightly smaller value of Za but which 
has the advantage of holding with an arbitrary magnetic field. It is this slightly weaker, 
but more verstile version of the stability of relativistic matter that enters the proof of the 
stability of the Pauli Hamiltonian with the field energy fffleid- (The ability to add a field 
was not explicitly stated in [13], but the extension is easy; it was first noted in [9].) The 
constants obtained in [13] have been improved recently in [14]. However, it remains true 
that the optimal value of Za — 2/7T has not been reached with a magnetic field present. 
The sufficient condition for stability in [14] (with a field) is 

(-)Z + 2.80 ^2/3 + 1.30 < 1/a , (3.3) 

which permits Z < 59 for a = 1/137. 

We are now ready to use (3.3) in the study of the stability of the nonrelativistic 
magnetic Hamiltonian (2.6). With Z fixed, let us define the number a by (3.3), i.e., 

1/a = (^)Z + 2.80 ^2/3 + 1.30 . (3.4) 

Then from (3.1) we conclude the operator inequality 

1 ^ 

V,>^y2\p + A\i. (3.5) 

This inequality can be substituted in (2.6), and we see that it now suffices to study the 
spectrum of the one-body operator 

a 

/iA = rA--|p + A| . (3.6) 

a 
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In fact, what we require is a lower bound on the the sum of the lowest N eigenvalues of 
/lA (according to the Pauli principle). One can show [9] that this sum is bounded below 
by const, x N — iffieid, as required for stability, provided that the conditions given in (2.7), 
or the ones mentioned after (2.7), are satisfied. 

An inequality that plays a key role in this proof is the CLR bound on the number of 
negative eigenvalues for a single particle in a potential (in n > 3-dimensions) , i.e., for the 
Hamiltonian 

h = p^ + U{x) (3.7) 

with = — as usual and with U some arbitrary potential. We can always write 
U in terms of its positive and negative parts, U{x) = U+{x) — U-{x), with C/+(a;) = 
{l/2){\U{x)\ + \U{x)} and U-{x) = {l/2){\U{x)\ - U{x)}. Then there is the inequality 
for the number of negative eigenvalues of h, obtained independently by Cwikel, Lieb and 
Rosenblum, 

Number of negative eigenvalues < Lo,n J U-{x)'^^'^ d^x (3.8) 

with Lo,3 = 0.1156. (The sharp value of Lo,3 is not known, but it is not much smaller than 
this value, obtained by Lieb.) 

Closely related to this is a family of similar eigenvalue inequalities, derived earlier by 
Lieb and Thirring [3], about the power sums of the negative eigenvalues, Sj of h, 

|e,-r < L^,n I C/-(a;)^+"/^ d'^x (3.9) 

e,<0 

for n > 1, all 7 > for n > 2 and all 7 > 1/2 for n = 1. This inequality, for n = 3 
and 7 = 1 is important in the proof [3] of the stability of nonrelativistic matter without 
magnetic fields. The case n = 3 and 7 = 1/2 plays a role in the study of -f^Dirac below. It 
is known that Li^g < 0.0404 and Li 3 < 0.06003. 



IV. Putting it all together: Relativistic matter with 
magnetic fields 

Up to now we have considered the nonrelativistic Schrodinger equation with a Zeeman 
term (i.e., with the Pauli operator) and showed that the magnetic field energy restores the 
otherwise lost stability provided ZoP" and a itself arc small enough (but well within the 
range of physical parameters). We have also considered the "relativistic" Schrodinger 
operator, with no Zeeman term, and found that stability requires that Za and a be small. 
The combination of the two might be expected to lead to insurmountable difficulties, and 
it was with some surprise that we found that these difficulties can be overcome if things 
are defined properly. 
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We shall try to combine the two by using a Dirac operator for the electron kinetic 
energy. Thus, we define the Dirac operator for a particle of mass m in a magnetic field by 

Pa = P< + 4 + /5™- (4.1) 

As usual, (p = a ■ O, where ctj and (3 are Dirac matrices. Our many-body Hamiltonian is 
formally similar to the one in (2.6), namely, 

N 

i^Dirac = J] (*) + OiV^ + i^fleld , • (4.2) 

1=1 

There is, however, one very important difference from (2.6): -ffDirac is unbounded 
below because the Dirac operator itself is unbounded below (recall that the Pauli operator, 
7a is not only bounded below, it is positive). According to Dirac the remedy is to fill all 
the negative energy states of the Dirac operator which (because of the Pauli principle) is 
the same thing as saying that the electron wave function must lie in the positive spectral 
subspace of the Dirac operator. 

There's the rub! Which Dirac operator are we talking about? There are at least two 
significant ones in our problem. One is the free Dirac operator T>q = + fim. The other is 
the Dirac operator with the magnetic field "Da- It would seem that the latter is the more 
important one since an electron can never get rid of its own magnetic field. Moreover, the 
choice X>0) is not gauge invariant, i.e., the multiplication of an electron wave function by a 
spatially varying phase usually takes a positive energy function into a mixture of positive 
and negative energy functions. The second choice is manifestly gauge invariant. (To avoid 
possible confusion it is to be noted that the Hamiltonian is always given by (4.2); the only 
question is what condition to impose on the electron wave function.) 

This issue is usually not clearly stated in field theory textbooks, but whatever one 
might think about the appropriateness of either definition the interesting fact is that 
"stability" of matter can settle the argument. The following two things are proved in 
[1]. They show that one choice is valid and the other is not. (Recall that a was defined in 
(3.4).) 

THEOREM: 

(i). If the electron wave functions are restricted to lie in the positive spectral subspace of 
the free Dirac operator then the Hamiltonian i^oirac is unstable. More precisely, for 
any given values of a > and Z > 0, there are (sufficiently large) particle numbers 
N and K for which the bottom of the spectrum of Hj^ire^c is —oo. 

(a). If the electron wave functions are restricted to lie in the positive spectral subspace of 
the Dirac operator Dp^ then i^Dirac is stable (i.e., i^Dirac > provided a and Z are 
small enough. In particular, if we define olc to be the unique solution to the equation 

1 - {ac/af = {SnLi^^aef/^ (4.3) 

then a < ac suffices for stability. 
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As a final remark, we draw attention to a recent result [15]. If we consider the Dirac 
hydrogenic atom without magnetic field, with Hamiltonian Vq — Za/r, and if the electron's 
wave function is required to lie in the positive spectral subspace of Vq, the critical value of 
Za, below which this Hamiltonian is bounded below, is raised from the value 2/7r (which 
is the corresponding value for the "relativistic" Hamiltonian in Sect. 3) to the value 
(ZQ;)crit — 2/(7r/2 + 2/7r). As yet it is unknown to what extent a magnetic field will change 
this value further. 
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